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Abstract 

We obtain the next-to-next-to-leading-log renormalization group improvement of 
the spectrum of Hydrogen-like atoms with massless fermions by using potential NRQED. 
These results can also be applied to the computation of the muonic Hydrogen spectrum 
where we are able to reproduce some known double logs at 0{mofi s ). We compare with 
other formalisms dealing with log resummation available in the literature. 
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In Ref. (see also 0), the renormalization group (RG) improvement of the Heavy 
Quarkonium spectrum for the equal mass case was obtained within the potential NRQCD 
(pNRQCD) formalism ||. This result was compared with the one of Ref. [Q (see also |5[ 
obtained within the vNRQCD formalism [|7|. Disagreement was found. This disagreement is 
potentially important as it propagates to different observables like, for example, t-i produc- 
tion near threshold, where it is claimed || that the resummation of logs plays an important 
role. For instance, the matching coefficient of the electromagnetic current, which is a neces- 
sary ingredient in these calculations, is different [f3, El. Nevertheless, for the known logs at 



next-to-next-to- leading |T0[ and next-to-next-to- next-to-leading order [TT[ , both calculations 
happen to agree with each other. 

In order to try to clarify this issue, we will consider the simplified problem of a Hydrogen- 
like system coupled to nf light (massless) fermions in QED. We will then obtain the next-to- 
next-to-leading-log (NNLL) RG scaling of the spectrum of this system. In principle, these 
results can be applied to muonic Hydrogen. In this case, the electron is replaced by the 
muon, nf — > 1 and the remaining light fermion is the electron (which we will take to be 
massless for simplicity or, at most, of 0(ma 2 ), where m is the mass of the muon). In this 
situation, we will be able to compare, in certain limits, with finite 0(ma 6 In 2 ) results already 
available in the literature [|12[]. Our results will agree with these calculations. 

The computation will closely follow the procedure of Ref. [1|] to which we refer for details. 
Here we will just write the main formulas necessary to set up the notation and the results. 

The first step is to obtain the RG improved matching coefficients of the NRQED ]13[ 
Lagrangian at one loop and up to 0(l/m 2 ) (m is the mass of the massive lepton (the muon 
for the muonic hydrogen) and the mass of the nucleus is sent to infinity in this paper). 

The NRQED Lagrangian including light fermions reads at 0(l/m 2 ) (up to field redefini- 
tions) 0, [TJ, |T§ 

£ = £ph + Ci + + C p + C^ p , (1) 

where fi is the Pauli spinor that annihilates the fermion, N p is the Pauli spinor that annihi- 
lates the nucleus, iD = id — gA , iD = i'V + gA, 



d = + E + c " j-2 E kfl^i T finish, (3) 
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ido + qZAq and 
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We have also included the D 4 /m 3 term above since it will be necessary in the evaluation of the 
spectrum once the power counting is established. Moreover, we will consider that the kinetic 
term matching coefficients are protected by reparameterization invariance (cj. = C4 = f) [16| , 
however, we will often keep them explicit for tracking purposes. 

By definition, NRQED has an ultraviolet cutoff z/nr = {u p , v s } satisfying mv <C z/nr <C m. 
u p is the ultraviolet (UV) cut-off of the relative three-momentum of the heavy fermion and 
antifermion. v s is the UV cut-off of the three-momentum of the photons and light fermions. 
The derivation of the scale dependence of the matching coefficients with respect the UV 
cutoffs of the theory goes identical to the one in Ref. [pi. In particular the fact that no 
dependence of v p appears at this order. In principle, the running of c and c M ' could be 
deduced from the results of [O, KM by taking care of the changes of the color structure. 



Since we are only interested in the computation of the spectrum at NNLL accuracy, their 
contribution will vanish at this order as far as the spectrum is concerned (c^ appears in the 
equation of cd but the running of is zero at LL accuracy). Therefore, the relevant RG 
equations in our case read 

d 



Vs-r-c D 



01 ( 8 r 2 , 00 al 

"tt \ 3 T 



and zero otherwise. 

By taking the matching conditions at the scale m: Ck 



(7) 



cf = c s = cd = 1 and {d} = 0, 



we can obtain the solution of the RG equations. We only explicitely display those which will 



be necessary later on (we define z 
T F = 1) 
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The above results are a necessary step towards the RG improvement of pNRQED with 
the matter content described above, which we consider in what follows. pNRQCD is defined 
by the cut-off z/ p nr = {v p , v us }, where v v is the cut-off of the relative three-momentum of the 
heavy fermions and is such that mv <^v v and v us is the cut-off of the three-momentum 
of the photons and light fermions with mv 2 v us mv. 



The pNRQED Lagrangian reads as follows (iDg = ido + g{Z 

'xd 3 X5 t (x,X,f) 
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where x and X, and p and P are the relative and center of mass coordinate and momentum 
respectively. All the gauge fields in Eq. are functions of the center-of-mass coordinate 
and the time t only. We have explicitly written only the terms relevant to the analysis at 
the NNLL. 

We now display the structure of the matching potentials V^°\ V^ 1 ' and V^ 2 \ which are 
the relevant ones to our analysis. At order l/m°, we have the static potential: 



y(0) ^ _ Z ^V_ (1Q) 

r 



In principle, at order 1/m, we may have a potential scaling as ~ — t-. Nevertheless, it 
vanishes at the order we are working. It would give, at most, 0(ma 6 ) corrections to the 
spectrum in a finite order calculation and the running equations would not mix with it. 
Therefore, for the purposes of this paper, we approximate 

~0. (11) 

m 

At order 1/m 2 , to the accuracy we aim at, has the structure 

= ^V (r) + 3^ii Li . (12) 

m 2 m 2 2m 2 r 3 

where Si = cr 1 /2. In principle, one may consider more structures for the 1/m 2 potential 
but, since they will not contribute at the accuracy we aim and in order to focus the problem 
as much as possible, we will set them to zero in what follows, as we have done for the 1/m 
potential. 

The coefficients, V = {aty s , D s , ...} contain some lnr dependence once higher order 
corrections to their leading (no n- vanishing) values are taken into account. In particular, 
we will have expressions like 5®(r) ln n r. This is not a well-defined distribution and should 
be understood as the Fourier transform of ln n 1/k. Nevertheless, in order to use the same 
notation for all the matching coefficients, and since it will be sufficient for the purposes of 
this paper, namely to resum the leading logs, we will use the expression d>( 3 '(r) ln n r, although 
it should always be understood in the sense given above. 

By studying the UV behavior of pNRQED it is possible to obtain the scale dependence 
of the coefficients of the potentials V. The discussion closely follows the one of Ref. |3J to 
which we refer for details. Here we just mention the main points. The potentials have the 
following structure: 

V(d{u p , u s , m), c(u s , m), u s , v us , r) = V(u p , m, u us , r) = V(v p , u us ) . (13) 

In particular, 

V »T-V = °- ( 14 ) 
dv s 

Moreover, at the accuracy we aim, we also get 

^ = °- (15) 
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Therefore, we obtain 

V(p p , v us ) ~ V{v us ) 

and we only have to compute the v us scale dependence. 

The i/ us -scale dependence could be obtained along the same lines as in Ref. 
obtain in this specific case: 

„2 
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and zero for the other potentials. 

Eqs. (P}), (|T3p and ([H]) provide the complete set of RG equations at the desired order. 
By using Eqs. ( |14"D and ([15|), we obtain 



V = V(d(l/r,m),c(l/r,m),i' s = l/r,v us ,r) . (18) 

We now need the initial condition in order to solve the US RG equations, i.e. the matching 
conditions. We fix the initial point at v us = 1/r. In summary, we need to know the 
static potential with 0(a 3 ) accuracy, the 1/m potential with 0(a 2 ) accuracy, the 1/m 2 
potentials with 0(a) accuracy and Va with 0(1) accuracy at v us = 1/r. For the non- 
vanishing potentials, they read 
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(19) 

where /3i = —ATpUf and the values of a x and 02 can be easily obtained from the QCD results 
by taking C f I, C A and T F -> 1. 

We now have all the necessary ingredients to solve the RG equations. The RG improved 
potentials read: 



otv\y us ] 



a v (r x ), 

dp 
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This completes the RG evaluation of the pNRQED Lagrangian at NNLL. 



(20) 
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With the above results we can obtain the energy with NNLL accuracy. The discussion 
goes similar to the one in Ref . [|lj . All the large logs can be obtained from the potential terms. 
Once the potentials are introduced in the Schrodinger equation, the ln ra (l/r) terms produce 
ln n (ma) terms plus subleading contributions (ln n_1 (ma), ■ • •) within the LL resummation 
counting. The expectation value of the potential terms is u us -scale dependent. This scale 
dependence is cancelled by the ultraviolet scale dependence of ultrasoft loops. The typical 
scales in these integrals is of the order ma 2 . Therefore, the logs of the ultrasoft loops 
get minimized by setting v us ~ ma 2 and all the large logs get encoded in the potential 
contributions. Finally, one obtains the following correction to the NNLO energy expression: 

^>J = ^(-|lo g (^)-3(c D -l)), (2!) 

where E n = —mZ 2 a 2 / (2n 2 ) and the scale v s in z and in the NRQED matching coefficients 
has been fixed to the soft scale v s = 2a" 1 , where a' 1 = mZa( - 2an \ a j s a i so understood at 



the soft scale v s = 2a" 1 unless the scale is specified. The z/ us -scale dependence of Eq. ( pl"D 
cancels against contributions from US energies. Since ma 2 is the next relevant scale, their 
effective role will be to replace v us by ma 2 (up to finite pieces that we are systematically 
neglecting) in Eq. (^T|). In particular, we take v us = —E n . As expected, Eq. ( pT]) with 
v us = —E n reproduces the well known Hydrogen-like 0(ma 5 In a) correction but, indeed, 
Eq. (pTD gives all the 0(ma 4 (a lna) n ) terms for n > 1 of the spectrum of the Hydrogen-like 
systems with rif massless fermions. After adding to Eq. ( pip the NNLO result with the 
normalization point at the same soft scale, v s = 2a~ 1 , that we have used here, the complete 
NNLL mass is obtained. Note that the above resummation of logs also correctly accounts 
for In Z terms with the same accuracy. 

We have seen that the large logs of the spectrum can be obtained from the potential 
terms by setting 1/r ~ ma and v us ~ ma 2 . The velocity of the non-relativistic particle is 
typically v ~ Za. Therefore, it is interesting to consider the scaling of the potentials with 
respect v as it will help us to later compare with vNRQED results. In practice, we will 
consider its scaling with respect v = mv (therefore v us = v 2 jm) where 

V{v p , m, v us , r) ~ V(m, u us ,r) -> V(m, v 2 /m, 1/V) = V(v) . (22) 

(2) 

We can now consider its derivative with respect v. We will just focus on D d since it is the 
only one which has a non-trivial running. We obtain 
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u—Df = -^ CD (v)a 2 (v) + -— ^ In -U- - --^a(- • 2 3 
du 47r 3 7r a(— ) 3 7r m 



■ m - 



It is remarkable that the above expression can be rearranged as 

d ^(2) Po ,v 2 s o, s 8a(v) ,v 2 . . 

u—Dy = -—c D (—)a 2 (v) —a(-) . 24 

dv d 4vr y m' v ' 3 tt y m' v ; 

There is an evaluation 0j within the vNRQCD framework of the RG improved Heavy 
Quarkonium mass when Aqcd ma 2 . The evaluation performed within the pNRQCD 
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framework [|Xj] disagreed with that evaluation. It was noticed there that the disagreement 
still persisted if one considered a QED-like limit with light fermions by taking Cf — > 1, 
Ca — > and Tp — > 1. Agreement was found for a QED-like limit without light fermions by 
taking Cf — > 1, Ca — > 0, n/ — > 0, Tp — > 1. Some errors seem to have been detected in the 
first versions of these calculations in vNRQCD (TB| that may partially explain the difference, 



in particular for the 1/m 2 potential. In this case, agreement may exist in the limit Cf — > 1, 
C A -»• and T F -> 1. 

For the evaluation performed in this paper, the computation of the spectrum for the 
case of Hydrogen-like atoms with massless fermions, there exists no analogous within the 
vNRQED framework. Nevertheless, it is possible to guess what would be the result in that 
formulation by using the rules of Ref . |19| , which relate the anomalous dimensions computed 



here with the ones that should appear in vNRQED. For the specific case of D d , we obtain 

v —D {2 ) (vNRQED) = 7s + 2 7n , (25) 

where 

Is = -—c D (v)a (u) , 7„ s = a(— ) . (26) 

47T 3 it m 

This should be compared with the running in pNRQED obtained above. If we do so, we 
find that Eqs. (p5j) and (24) are different. If expanded in a they first differ at 0(a 2 log 2 a). 



This produces a difference in the computation of the mass at O (ma 6 log 2 a). In order to 
perform an independent check, it would be extremally important that corrections of this 
order had been computed before. The closest system to the one discussed here corresponds 
to the muonic hydrogen for which, indeed, corrections to the energy at this order have been 



computed by Pachucki [12]. In order to compare our results with his evaluation, we have to 
take the limit n/ — ► 1. Moreover, for the real muonic hydrogen, the mass of the light fermion 
(the electron in this case) is not negligible. However, we can formally consider the situation 
m e ~ ma 2 (even if for the physical situation m e ~ ma is closer to reality) in his and our 
calculation. For the matter of comparison, in our case, this means that, for scales of the 
order of m e and ma 2 , we can use the low energy electromagnetic coupling a em ~ 1/137.... 
This is indeed the parameter expansion used in Pachucki's calculation. A closer inspection 
shows that the diagrams that give rise to the large logs computed here correspond to the 
ones drawn in Fig. 4 in Ref. |I2" |. If we reexpand our result in terms of a em = a(i/ MS ), we 



obtain (up to the order of interest and with v ~ a) 

n ( 2 ) a{u us ) f /3 . mv 2 \ 

° d ~— = —{ 1+ 2n- aMln ^ + ---) (27) 

8a(u us ) mv 2 2 fa(u us )\ 2 2 mv 2 

In -A) In + ■■ 

6 7i m 3 y 7r j m 

4a 2 (i/« s ), m v 2 2/3n a 3 (zA, s ) , mv 2 , mv 2 1 „ a 3 (^ ns ) , o mv 2 

~ ^^ln — — V^ln In f3 — In 2 . 

3 7T m 3 7r 2 mv m 3 n 2 m 

It is easy to identify the above terms (last equality) within a diagrammatic picture. The 
first term is the standard Lamb-shift correction one would find for the Hydrogen atom and 
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corresponds to the diagrams of Fig. 4 of Ref. JT2J without any bubble insertion. The second 
term corresponds to the first diagram in Fig. 4 of Ref. [12j[. The last term corresponds 



to the second diagram in Fig. 4 of Ref. 



Therefore, our result seems to have the 



correct structure for the 0(ma 6 ln ) corrections. Let us now go deeper in the comparison 
with Pachucki's results. First, we can see that the last term of Eq. fl27p can reproduce the 
analogous Pachucki's contribution by setting a{v us ) = a em and v us = m e (this result depends 



on the two-loop muon form factor first computed in Ref. [2D|). For the second term of Eq. 
(PT|), the explicit comparison is a little bit more involved. Nevertheless, it is possible to see 
that the first term in Eq. (39) of Ref. 



12] gives the logs of the second term in Eq. 



one can make the replacement (as far as the LL contribution is concerned) 



V VP 



Za(v v 



A) / mv 2 

—a{u us ) In 

2tt mv 



since 



(28) 



for VyPi as defined in Ref. 0§. The second term in Eq. (39) of Ref. |12j gives the logs due 
to expanding the wave-function at the origen ~ (ma) 3 (which are naturally written in terms 
of ot(v)) in terms of a em . Therefore, we can trace back all the logs of the computation in Ref. 
|i~2|l . This provides a check of our calculation to a level where it starts to first differ with 
what would be the vNRQED result. Nevertheless, it may happen that, if the corrections 
of the vNRQCD results for the equal mass calculation are finally confirmed, they may also 
explain the different result obtained here. 

In conclusion, we have computed the energy spectrum at NNLL for an Hydrogen-like 
system with rif massless fermions. We have checked our results at 0(ma 6 In 2 ) by comparing 
with results already available in the literature [12] for muonic Hydrogen and agreement has 
been found. We have also compared with what we would expect to be the result in the 
vNRQED framework based on the rules of [I9~ and disagreement has been found. Finally, 
we would like to mention that the above results can be useful in checking higher order logs in 
computations of the spectrum for muonic atoms or alike where the electron can be considered 
to be a light particle. 
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Note added. After this paper was sent to the hep-ph archives a new paper [hep-ph / 0204299 
appeared in the net where it was pointed out that there was a systematic error in the 
vNRQCD computations to date and that the diagram Fig. 20b in this reference should be 
included in such computations. After its inclusion, the corrected 1/m 2 potential obtained 
in vNRQCD agrees with the 1/m 2 potential obtained in pNRQCD PJ, the replacement 
Cd{v) — > Cr){v 2 /m) should be done in Eq. ( p6|) and the new result of vNRQED for the 
muonic hydrogen spectrum agrees with our result. 
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